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Chapitre 3 

Equations aux différences 

3.1 Equations aux différences 

3.1.1 Introduction et motivations 

Plusieurs algorithmes d'analyse numérique (résolution numérique d'équations différentielles) nécessitent 
de trouver les solutions d'équations aux différences. De plus, les équations aux différences jouent un 
rôle important dans d'autres branches de mathématiques pures et appliquées telles que: recherche 
des zéros des polynômes, théorie des probabilité et les mathématiques économiques. 

3.1.2 Equations aux différences du premier ordre 

Définition: 

Soit {f n (y, z)) n çjN une suite de fonctions, y, z G S C 1R. Le problème est de trouver une suite (x n ) 
telle que: 

(i) x n G S , X n -l G S 

(ii) f n (x n ,x n -i) = équation aux différences 

Une suite satisfaisant (i) et (ii) est appelée solution de l'équation aux différences (ii). 

Exemples: 

a. La suite f n (y,z) = y — z — 1 donne l'équation aux différences x n — x n -\ = 1. La solution de 
cette équation s'obtient en écrivant: 

%n— l X n —2 — J- 

x l ~ x — 1 
en faisant donc la somme on trouve: x n = n + xq. 

b. La suite f n (y,z) = y — z — n, donne l'équation aux différences x n — x n _i = n. La même 
technique que celle utilisée en a.) donne: 

n(n + 1) 

x n = " h X 

35 
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36 2. Equations aux différences 

3.1.3 Equations aux différences d'ordres N 

Définition 

L'équation aux différences d'ordre N (N G IN*) est définie de façon similaire. Soit {fn(yo-, yi-, • • • •> Vn)) 
une suite de fonctions avec g/Q,2/i, . . - , ijn é^cE. Le problème est de trouver (x n ) telle que: 

(i) x n € S , x n _i € S > . . . , x n -u 6 S 

{iï) f n (x n ,x n _i, . . . ,x n _j\[) = équations aux différences 

Une suite satisfaisant (?) et (il) est appelée solution de l'équation aux différences d'ordre N. 

Exemple 

La suite de fonctions définie par: f n (y, z,t) = y — z — yt donne l'équation aux différences suivante: 
x n ~ x n-i ~ x n-i x n-2 = 0. C'est une équation aux différences non linéaire. 

3.2 Résolution des équations aux différences 

On se limite au cas linéaire. 

3.2.1 Equation du premier ordre 

Soit à résoudre x n = a n x n -\ + b n avec xq = cq = c. 

• première étape: équation sans second membre: x n = a n x n -\. 
En écrivant cette équation pour n. n — 1 . . . , 1 on obtient: 

Xn — QinXji— 1 
X n —\ — Q, n —\X n — i 

X\ = a\XQ , 
en multipliant terme à terme on obtient: 

X n = O1O2 ■ ■ ■ fln-lfln/ = n n ;r n 

n„ 
On pose x n = cli n , Llq = 1 

• deuxième étape: équation avec second membre. 

On utilisera pour celle-là la méthode de la variation de la constante: 

x n = c n Il n , où c n est à déterminer 

En injectant x n = c n W n dans l'équation aux différences avec second membre x n = a n X n -\ + b n on 
trouve: 

C n Ii n = a n Cn-lKn-l + b n = Cn-lUn + K 

On suppose que a n ^ (n = 1,2,... ), donc n n 7^ 0, la suite (c n ) vérifie donc l'équation aux 

différences suivante: 

, b n 

Cn = c n -i + — - 
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2. 2 Résolution des équations aux différences 37 

en écrivant cette équation pour n, n — 1, . . . , 1, on obtient: 

, h. 

Un. 

bn-1 

Cn-1 = C„_2 + :pj 

■II71-I 

Cl = C(, + nT' 



En faisant la somme on trouve: 






la solution générale est donc: 



Théorème: Soit a n 7^ 0, n = 1,2,.. .La solution de l'équation aux différences x n = a n x n _i + b n 
telle que xq = c est donnée par: 

X n = ïln{C + — + ...+ - + =-) 

lll H71-I Ain 

où ITo = 1 , II n = aia 2 . . . a n 

3.2.2 Equations aux différences du second ordre 

Les équations aux différences du second ordre sont de la forme: 

X'n T O n X n — 1 -p C n X n —2 — "n i "ras '-•ri: ^*ti t vb 

nous donnons la solution générale de cette équation aux différences dans le cas où d n = 0, b n = b 7^ 

et c n = c 7^ 0. 

Soit £ l'ensemble des suites de nombres réels vérifiant: 

x n + bx n _i + cx. n _ 2 = . (3.1) 

On montre facilement que £ est un sous espace vectoriel de l'ensemble des suites réels de dimension 
2 sur H. 

On se propose de chercher les solutions de la forme x n = r n où r est un réel. Si x n = r n est 
solution de 3.1 alors r vérifie: 

r 2 + br + c = (3.2) 

Trois cas se présentent: 

• A = b 2 — 4c > 0, dans ce cas 3.2 admet deux racines réelles distinctes ri et T2- La famille (r™ , r^) 

est une base de £. La solution générale de 3.1 est donnée par: 

x n = ar^ + 0r% 

a et (3 peuvent être calculés facilement en fonction de xq et x\ et sont donnés par: 

$1 ~ XqV 2 ,, Xon - X\ 

a = , p = 
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38 2. Equations aux différences 

finalement on trouve: 

_ xi - x Q r 2 n xpn - xi n 

x n — ri -\ r 2 

ri — r 2 r\ — r 2 

• A = 0, l'équation 3.2 admet une racine double r^. Dans ce cas les deux suites r^ et nr^ forment 
une base de £. La solution générale de 3.1 est donc: 

x n = x r^ + {z\ - x {) r d )nry l 

• A < 0. dans ce cas 3.2 admet deux racines complexes Ai = pe et Ai = pe . On montre 
également que les suites p n cos(n9) et p n sm.(nô) forment une base de £. La solution générale est 
donnée par: 

f û\ , Xl ~ x 0P cosB n ■ i m 

x n — ^oP coslnv) H p sin ra 

p sin 9 

3.3 Schéma de Horner 

3.3.1 Algorithme 1 

Soit à calculer xq, Xi, ... xn de façon récursive à partir de: 

Xq = bg , x n = zx n _i -\-b n , n = 1, 2, . . . , N 

c'est une équation aux différences du premier ordre dans un cas spécial: a n = z et donc II n = 
z n ,n = 1,2, ...N. Sa solution est donc donnée par: 

x n = U n (b + — + ... + -^-) 

tll i-i-n 



Théorème 

Les nombres produits par l'algorithme précédent sont égaux à la valeur du polynôme P n défini par: 

P n (X) = b X n + hX 11 ' 1 -f . . . + b n ,n = 0, 1, . . . , N 

pris h X = z. 

L'algorithme précédent est appelé schéma d'Horner. 

3.3.2 Nombre d'opérations: d'additions et multiplications: 

Soit P n un polynôme de degré n. 

Quelle est le nombre d'additions et de multiplications nécessaires à effectuer pour évaluer P n (z)7 

P n {z) = b Q z n + hz 11 ' 1 + . . . + bn-xz + b n 

• évaluation directe: 

Chaque terme bz l nécessite i opérations de multiplications, P n (z) nécessite n opérations d'additions. 

Le nombre total est donc: 

/V^ -x n ( n + 1) 



évaluation en utilisant 



~n ~n— 1 , 
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2. 3 Schéma de Horner 39 

A l'aide de cette relation on remarque que si z 77-1 est connu, alors z n nécessite 2 opérations de 

multiplications. En écrivant: 

P n (z) = b n + b n _iz + b n _ 2 zz + b n -zz z + . . . + b x z n ~ 2 z + b G z n ~ 1 z 

L'évaluation de P n {z) nécessite donc 2n— 1 opérations de multiplications et n opérations d'additions. 
Soit au total 3/i — 1 opérations. 

• évaluation en utilisant le schéma d'Horner 

L'évaluation de P n (z) à l'aide du schéma d'Horner est donnée par l'algorithme suivant: 

x = b 

x\ = b\ H- zxq 

x 2 = b 2 + zx\ 

Xyi n ~\- ZXji—1 

ceci est équivalent à: 

Pn(z) = b n + z(b n -i + z(b n _ 2 + • • ■ + z(bi + zb ))) 
chose que l'on peut schématiser comme suit: 



le triangle 



signifie que Xi = bj + zXi—\. Il est facile de se convaincre que x^ = Pn( z )- 

D'après l'algorithme, on remarque qu'il faut n additions et n multiplications pour évaluer P n (z), 
soit au total 2n opérations. 

Conclusion: 2n < 3rc — 1 < n{n + l)/2 + n donc le schéma le plus rapide est le schéma d'Horner. 

Exemple 

Soit à calculer P(X) = 7X 4 + bX 3 - 2X 2 + 8 en x = 1/2. 

On a 6 = 7 , h = 5 , b 2 = -2 , b 3 = , 64 = 8 et donc: 

Xi = 8.5 , x 2 = 2.25 , x 3 = 1.125 , .t 4 = 8.5625 

3.3.3 Application du schéma d'Horner à l'évaluation de la dérivée d'un polynôme 

Soit P(X) = b X N + biX 1 *' 1 + . . . + b i X N ~ i + ... + b N un polynôme de degré N. Le développement 
de Taylor de P au voisinage de X = z s'écrit comme: 

P(X) = P(Z + h) = C h N + dh 1 *- 1 + . . . + Cih 1 *- 4 + . . . + Cyv 

avec h = X — z et 

CN-k = ^P (k) (z) ; k = 0,l : ...,N 
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40 2. Equations aux différences 

3.3.4 Algorithme 2 

Soit (xn ) la suite générée par l'algorithme 1 pour k = 1,2, . . . , N. Nous définissons la suite (xn 
récursivement par: 

4 t) =4 fc - 1) , x« = ^2 1 + 4*- 1 ' , n = 0,l,2,...,N-k 

k = , n = 0,1,2, ,N 

k = , n = 0,1,2, ,JV-1 

fc = JV - 1 , n = 0, 1 

k = N , n = 

Exemple 

N=4 

x (o) x (o) x (o) x (o) x (o) 

I 

x (1) 

Xq 

I 



x (3) 

Xq 
r W 

Xq 

avec 

4 0) = 6oz 4 + M 3 + hz 2 + & 3 z + h 

Théorème 

Soit P n la suite des polynômes: 

P„{X) = b a X n +b 1 X"- 1 + ... + &„ , n = 0,l...,JV 
Les suites (a4 ) de l'algorithme 2 vérilient: 

r' fc ' = — P^' r, = 12 AT £- = 01 AT 

x„ — ^ r n+ k i" J-i^i • • • >-" i n. — u, i, . . . ,iv 

Si -P(^) = -P(-z + /l) = Co/i + c\h + . . . + cat les coefficients CN_k sont les éléments diagonaux 
du schéma suivant: 

r (<» r (0) r (o) (0) (o) (o) 

(1) (1) (1) (1) (1) 
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4 N) 



^ETtfUP 



2. 4 Application à la recherche des zéros d'un polynôme 
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Preuve: par récurrence sur k (très technique et est laissée à titre d'exercice pour les étudiants 
ambitieux [Exercice 8]). 

Exemple 

Ecrire le développement de Taylor de P(X) = 7X 4 + 5X 3 - 2X 2 + OX + 8 au point x = 4: on a 
b = 7 6i = 5 b 2 = -2 b 3 = b i = 8 
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3.4 Application à la recherche des zéros d'un polynôme 

3.4.1 Méthode de Newton— Raphson 

Nous allons exposer cela à travers l'exemple suivant: Soit P(x) = x 3 — x 2 + 2x + 5 un polynôme 
de degré 3. On a P{— 1) = 1 et P(— 2) = —11, donc le polynôme P admet une racine s G [—2, — 1], 
Nous cherchons une approximation de s en utilisant la méthode de Newton— Raphson: 

xq donnée 

T 1 -T - P ^") 
■^n-r-l — -^îl P'(x ) 

à partir de xq, nous calculons les autres termes de la suites (x. n ). P(x n ) et P'(x n ) seront évalués 
en utilisant le schéma d'Horner, on a donc le tableau suivant: 



n 


Xfi 


P(x n ) 


P'(x n ) 





-1.5 


-3.625 


11.75 


1 


-1.191489 


-0.4941198 


8.641919 


2 


-1.1343122 


-0.01476801 


8.12861753 


3 


-1.13249547 


-0.00001452 


8.1126289 


4 


-1.1324936884781 


-1.41 10~ n 


8.112613 


5 


-1.1324936884764 







3.4.2 Méthode de Bernoulli (1654-1705) 

La méthode de Bernoulli permet d'obtenir la racine de plus grand module d'un polynôme. Pour 
cela, on associe à l'équation polynômiale une équation aux différences. Le rapport de deux ter- 
mes successifs d'une suite, solution particulière de l'équation aux différences, tend vers la racine 
cherchée. Avant d'ennoncer le théorème général, nous étudions l'exemple suivant: 



Exemple: Soit l'équation x - 

différences associée: 



5a: + 6 = dont les racines sont 2 et 3. Considérons l'équation aux 



x n+2 ~ 5x n+ i + Qx n 



(3.3) 
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42 2. Equations aux différences 

Comme nous l'avons montré à la section 2.2.2, les suites u n = 2 n et v n = 3" sont solutions 
particulières de 3.3. La solution générale est donnée par: 

x n = a2 n + (3T , a,^6H (3.4) 

On peut toujours choisir Xq et ^i de telle sorte que j3 ^ 0. 
Considérons le rapport Xn+1 , d'après 3.4 ce rapport s'écrit comme suit: 

x n+1 Q 2" +1 + /33" +1 /33" +1 (f (|)" +1 + 1) 



n2" 



-/33" /33"(|(|)" + 1) 



On remarque que: 



(f CD- + 1) ( 5) 



lim ^±1 = 3 



— = E 




il 




ïî - 


3.8 


X 2 




X4 


3.42105 


%3 




x 5 


= 3.24615 


X4 




x 6 


= 3.15165 


Xç 




x 7 


= 3.09624 


X(, 




Z8 


_ q n«oi fii 



Application numérique: 

Si nous choisissons comme conditions initiales: xq = 0, x\ = 1 alors 

x 2 = 5 

2:3 = 19 

2:4 = 65 

2:5 = 211 

xg = 665 

x 7 = 2059 

2 S = 6305 , — = 3.062166 

x 7 

On constate sur cet exemple que le rapport x " +l tend vers 3 lorsque n tend vers l'infini. 

3.4.3 Théorème général: 

Soient m un entier positif donné et (ûi)i=o,....m rn + 1 nombres complexes donnés. On considère 
l'équation aux différences: 

O-O^n+m + a-lXn+m-1 + a,2X„ +nl -2 + ■ • • + a m x n = (3.6) 

i) L'ensemble des solutions de l'équation 3.6 forme un espace vectoriel de dimension m 

ii) Si les racines ri, r 2: . . . r m de l'équation 

a x m + ai 

sont simples, alors la solution générale de l'équation aux différences 3.6 est de la forme: 

x n = a ir 1 + a 2 r^ + a-^ + ... + a m r™ 



^ETtfUP 



2. 4 Application à la recherche des zéros d'un polynôme 43 

iii) Soit P(x) = aox m + a\X m + . . . + a m un polynôme à coefficients complexes. Supposons que 

toutes les racines (r^ = i m sont distinctes et qu'il existe une racine ri telle que |n| > |r$| 

î = 2, 3, . . . , m. Alors si nous fixons xq, Xi, . . . , x m —i de manière a ce que ai ^ 0, la suite 
(x n ) solution de 3.6 vérifie: 

-,. x n +\ 

lim = ri 

iv) La méthode de Bernoulli, dans le cas de racine dominante, est d'ordre 1. 
Preuve: 

i. Soient (x n ) et (y n ) deux solutions de 3.6, et soient a et /? deux nombres complexes, alors la 
suite (^J définie par: z n = ax n + /3y n est aussi solution de 3.6. Comme nous l'avons vu dans 
les exemples précédents, la solution (x n ) d'une équation aux différences du premier ordre (resp du 
second ordre) est complètement déterminée par la donnée de xq (resp xq et x±). De même, une 
solution particulière de 3.6 est entièrement déterminée par la donnée de Xq, X\, . . . , x n —\. 

ii. De façon similaire à l'équation aux différences du second ordre, on cherche des solutions de 
3.6 de la forme {r n ). Si c'est ainsi, alors r vérifie l'équation suivante: 

a r n+m + air" 4 ™" 1 + a 2 r n+m - 2 + . . . + a m r n = (3.7) 

a (r n r + a^r- 1 + a 2 (r") m " 2 + . . . + a m r n = (3.8) 

comme le montre l'équation 3.8, les solutions non nulles (de la forme r n ) sont les racines de l'équation 
polynômiale: 

a () x 7n + ai x m ~ l + . . . + a m = (3.9) 

Si les racines de 3.9 sont simples, alors, et du fait que l'ensemble des solution de l'équation 3.6 est 

un espace vectoriel, la solution générale de 3.6 s'écrit: 

x n = air" + a 2 r 7 2 + a z r^ + . . . + a m r™ 
avec ai e(C 
iii. A l'équation polynômiale P(x) = 0, nous associons l'équation aux différences: 

aoXn+m + aiXn+m-l + «2^n+m-2 + ■ ■ ■ + 0, m X n = 

comme nous l'avons signalé dans i, la donnée de Xq, X\, . . . , X m -\ permet de déterminer entièrement 
la solution (x n ) de l'équation aux différences précédente. On peut toujours choisir Xq, Xi, . . . , X m ~i 
de manière à ce que la solution x n = cx.\r^-\- 0.2^2 + Q 3 r 3 +■ ■ -+CtWm so ^ telle que ai 7^ 0. Calculons 
le rapport " T " +i : 

x n+ i «ir;' +1 + q 2 r 2 " +1 + q 3 rg +1 + . . . + Q»,^ 1 



X„ Qirf + 2 ' , 2 l + «3^ + ■ • • + "m^ 

n+1 1 J_ a ? I r ? \™+l _L Q 3 I r :i \n+\ 1 1 a m ( T m \n+l 

1 i ' apri' ûi Vri / ai v ri ' 

d'après les hypothèses émises, pour i ^ 1 on a : |^| < 1 et lim^-^+oo | — | n = 0, on conclut donc, 
d'après l'égalité 3.10, que: 

lim = ri 
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44 2. Equations aux différences 



iv. on se place dans les notations de iii, calculons l'erreur e n = Xn+i — r±. En utilisant la même 
technique que celle exposée précédement on trouve: 

Q2l>2 - n)rg + a[i(r 3 - ri)r% + . . . + QL m {r m - r{)r^ 
air? + «2^2 + a 3^3 + ■ ■ ■ + a m r^ 

= *(!ï)»l±^ (3.11) 

avec /% = oti(ri — r{) % = 2, . . . , m et 

P3 / r 3 v fl . . Pm / r rn v n 

^, = £«(!»)» + ^ ( îï r + ... + £!!» (!«•)». 

ai ri ai ri ai r\ 

Supposons que l'on a: \n\ > |ï™2 J > . . . > |r m | alors lim n ^ +00 ^ n = liirin— »-)- 00 ip n = 0. On a donc: 
1 e»+l 1 r 2 / (1 + ^n+i)(l + V J n) 



lim ,£ii±i| = T» i im f + ^ + W ) 



ra— ►+00 



ri 



donc la méthode de Bernoulli est d'ordre 1. 
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3.5 Exercices 

Exercice 1. Soit —1 < a < 1. Résoudre l'équation aux différences: 

x = 1 , x n = ax n -i + 1 
Déterminer la limite de la suite x n lorsque n — > +oo. 

Exercice 2. Soit z un nombre réel non nul. Résoudre l'équation aux différences 

n 
x = 1 , x n = —x n -i + 1 
z 

Montrer que: 

lim — -x„ = e z 

Exercice 3. Soit la suite aux différences du premier ordre définie par: 

x = 1 , x n = {yJn)Z~ 2 x n _i + — avec Z G H 

(277.)! 

a. Résoudre cette équation aux différences. 

b. Quelle est la limite l(Z) de ^Lx n lorsque n — » +oo. Donner une valeur approchée de 1(1) 

Exercice 4. Soit l'équation aux différences du premier ordre définie par: 

1 1 ^ 

^o — 1 -, x n — ~ x n-i H r avec Z ÉE H 

z n\ 

a. Résoudre cette équation aux différences. 

b. Quelle est la limite l(z) de z n x n lorsque n — > +oo? 

c. Donner une valeur approchée de 1(1) à 10 . 

Exercice 5. Supposons que la série X^J^Lq ^n es "t convergente, et soit s sa limite. Montrer que 
s = lim^^oo x n ou (x n ) est la solution de l'équation aux différences suivante: 

x = b , x n = X n -\ + b n 

Exercice 6. En utilisant le schéma d'Horner. trouver le développement de Taylor de P(x) = 
Ax 3 - 5x + 1 

Exercice 7. Soit P un polynôme de degrés n > 1: 

P(z) = a z n + a x z n ~ l + ... + ... + a n _iz + a n 
Soit Z\ une racine de P, on a donc la décomposition suivante: 

P(z) = (z - z 1 )P 1 (z) 
Montrer que Pi peut se mettre sous la forme: 

Pi (z) = boz"- 1 + b x z n ~ 2 + ... + ... + b n _- 2 z + b n -i 

ou les (bi) vérifient une équation aux différences que l'on précisera. 
Application: P(z) = z 3 - z 2 + 2z + 5 déterminer P\(z) pour z\ = -1.132494 
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46 2. Equations aux différences 

Exercice 8. Soit P un polynôme de degré n > 1; 

P(z) = a Q z n + a x z n ~ l + ... + ... + a n _iz + a n 
a.) Montrer que pour tout k = 0, 1, . . . , n on a: 

P<*>( z ) = kl(Cja ;? B -* + C^_ 1 a 1 z"- k - 1 + ... + ... + C^an-k^z + C h k a„_ k ) 
b.) Montrer que : 

fc "+" ^fc+1 "T ■ ■ ■ T O n — O n+1 

c.) En utilisant a.) et b.) montrer que la suite a4 définie par: 

r (*0 _ r (*-l) T (fc) _ z Jk) , (fc-l) 

Ou Xn est définie par l'algorithme 1, k = 1,2, ... N et n = 0, 1, 2, . . . , N — k est donnée par: 

X n - k ^n+k\ Z ) 

Applicationl: Calculer P(z) t P'{z) et P"{z) à z = 1.5 pour P(z) = 2z 5 - 7^ 4 + 3^ 3 -6z 2 + 4^- 5 
Application2: En ut lisant l'algorithme de Newton-Raphson déterminer la racine positive du 

polynôme P(z) = z 3 — z 2 — z — 1 
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